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When  a  liquid  jet  impacts  a  pool  containing  the  same  liquid  and 
surrounded  by  a  still  gas,  a  surface  depression  is  produced.  The  surface  shape  is 

1  2 


determined  by  the  Weber  number,  We  = 


2  Plvt 


# 

In  this  work  the  shape  of  the 


surface  is  obtained  as  a  function  of  the  Weber  number  by  using  an  asymptotic 
expansion  technique. 


INTRODUCTION 

The  entrainment  of  non-condensible  gases  by  a  plunging  liquid  jet 
impacting  a  liquid  pool  is  related  to  some  important  practical  problems.  For 
example,  the  absorption  of  greenhouse  gases  into  the  ocean  has  been  hypothesized 
to  be  highly  dependent  upon  the  air  carryunder  that  occurs  during  breaking 
waves,  which  represent  a  type  of  plunging  jet  [Monahan,  1991;  Kerman,  1984]. 
Other  applications  include  some  type  of  liquid/gas  chemical  reactors.  In  order  to 
enhance  the  reaction  rate,  a  jet  of  liquid  entrains  the  surrounding  gas  reactant 
forming  a  two-phase  jet.  The  reaction  rate  is  enhanced  because  of  the  increase  of 
the  interfacial  area  density  and  the  pseudo-turbulence  produced  by  the  entrained 
gas  bubbles. 

Other  less  obvious  applications  are  associated  with  the  slug  or  plug  flow 
regimes.  For  example,  in  vertical  slug  flow,  the  Taylor  bubbles  move  upward 
with  an  almost  constant  velocity.  The  liquid  film  surrounding  the  Taylor  bubbles 
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drains  downward  forming  an  annular  liquid  jet.  Therefore,  a  gas  entrainment 
problem,  similar  to  the  one  studied  in  this  work,  occurs  behind  the  Taylor  bubbles. 
In  particular,  a  Helmholtz-Kelvin  instability  occurs  at  the  interface  in  the  rear 
portion  of  the  Taylor  bubble  such  that  smaller  gas  bubbles  are  entrapped  in  the 
liquid  plug  behind  the  Taylor  bubble.  A  difference  with  respect  to  this  work  is  that 
in  the  slug  flow  case  the  velocity  distribution  is  nonuniform  while  in  this  work  a 
uniform  velocity  distribution  was  assumed. 

DISCUSSION 

Most  prior  theoretical  studies  were  done  for  liquid  jets  with  very  low  liquid 
velocity,  which  had  applications  for  the  coating  of  fiber,  etc.  Recently,  Lezzi  & 
Prosperetti  [1991]  performed  a  stability  analysis  considering  the  liquid  jet  and  the 
pool  to  be  inviscid  liquids  and  the  gas  to  be  viscous.  They  did  not  analyze  the 
surface  depression. 

Bonetto,  et  al  [1992]  analyzed  an  inclined  plunging  liquid  jet.  They  assumed 
that  both  the  gas  and  the  liquid  were  inviscid  and  they  obtained  the  entrained  gas 
flow  rate.  The  entrained  gas  volumetric  flow  rate  is  one  of  the  most  important 
quantities  which  we  want  to  compute  in  problems  of  this  type.  The  key  parameter 
in  such  evaluations  is  a  knowledge  of  the  gas  gap  thickness  (ie,  the  shape  of  the 
inclined  surface  depression. 

The  objective  of  this  work  was  to  evaluate  the  induced  surface  depression 
caused  by  an  impacting  liquid  jet.  We  have  assumed  that  both  fluids  are  inviscid 
and  irrotational,  and  that  the  gas  region  is  at  constant  pressure.  Hence  the 
appropriate  equation  for  this  problem  is  the  Laplace  equation  for  the  liquid  and 
the  associated  interfacial  jump  condition.  The  former  equation  has  no  parameter 
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and  the  latter  has  only  the  Weber  number,  We  = - — - ,  where  x0  is  the  half¬ 

width  of  the  liquid  jet. 

We  may  use  an  asymptotic  expansion  technique  to  solve  the  problem.  That 
is,  we  may  expand  the  solutions  in  terms  of  relatively  small  We  number,  and 
substitute  these  expressions  into  the  Laplace  equation  and  the  interfacial  jump 
condition. 

Equating  terms  of  the  same  order  in  We,  one  obtains  a  recursive  system  of 
equations,  which  can  be  numerically  evaluated  up  to  third  order.  In  other  words, 
we  can  evaluate  the  eigenfunctions  of  the  problem.  Once  we  have  these 
eigenfunctions  we  can  compute  the  shape  of  the  surface  (and  the  velocity  fields)  as 
a  function  of  the  We  numbers. 


ASYMPTOTIC  EXPANSION 

The  purpose  of  this  section  is  to  compute  the  position  of  the  interface,  q(x) 
(shown  schematically  in  Figure- 1). 

For  the  assumption  of  irrotational,  invisdd  flow  the  governing  equations  for 
the  liquid  field  are 

V2  v  =  0  (1) 


where  y  is  the  stream  function.  The  two  velocity  components  are  evaluated  as, 
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where  u  and  v  are  the  velocity  components  along  thex  and  y  axis,  respectively.  If 

A  Jvf 

we  prescribe  a  value  of  y,  or  its  derivatives  normal  to  the  surface,  :rr,  for  every 

ou 

point  on  the  boundary  we  have  a  well  posed  problem.  As  shown  in  Fig.-2,  the 
plane  x  =  0  is  a  symmetry  plane  then  the  lateral  u  velocity  must  vanish  for  every y 
atx  =  0.  Using  Eq.  (2) 

u(x  =  0,y)  =  |^(x  =  0,y)  =  0  (4) 

on  the  centerline  of  the  plunging  liquid  jet.  We  know  that  the  v-velodty  must  bev^ 
at  positions  where  the  jet  is  impacting.  From  Figure-1  and  Eq.  (3),  we  see  that: 

v  (0  <x  <  x0,y  *s  0)  =  -|^  (0  <x  <  x0,y  =  0)  =vt  (5) 

For  x  >  x0,  the  free  surface  position,  fj(x),  must  be  coincident  with  a  streamline. 

Without  loosening  generality  we  make  the  free  surface  coincident  with  the 
streamline  y  =  0,  then 

y(x,fj(i))=0  (6) 


We  may  assume  that  the  pressure  in  the  gas  region  is  constant  and  equal  to  zero. 
The  liquid  pressure  right  under  the  surface  is  related  to  the  curvature  of  the 
surface  by: 


°  dx2 

.  ^vy2 = P<  (x  .n  (x))  =  -  2  Pi  [“28  KM  +  v2£A  (*))] 


i+f^T 

Idx  I 


(7) 
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where,  a,  P^,  P;,  q  are  the  surface  tension,  liquid  pressure,  liquid  density  and 
position  of  the  interface,  respectively,  andu2(x,q(x))  and  v2(x,q(x))  are  the  u  and  v 
velocity  components. 

In  order  to  have  a  well-posed  problem  we  may  specify  y  for  £  — »  «>  for  every 
y,  V  (£  — >  «>y)  and  V  fs,y  — *  °°).  One  possible  set  of  boundary  conditions  is: 


V  f£  ooy)  =  o 

(8) 

V  (£,y  ->°°)  =  0 

(9) 

We  make  Eqs.  (l)-(9)  nondimensional,  using  x0,  v^  as  the  length  scale  and  velocity 

scale,  respectively.  It  is  convenient  to  work  with  nondimensional  quantities. 
Thus  we  have: 

v2V  =  0  (10) 


with  boundary  conditions, 


(x  =  0,  y)  =  0 

(11a) 

V  (x,  y-^«>)  =  0 

(lib) 

V  (x— >oo,y)  =  0 

(11c) 

|*<x,y=0)  =  -l 

0<x<  1 

(lid) 

V(x,q  (x))  =  0 

1  <x 

(He) 

where, 


u  =  x/x0 


y=y/*0  *  ii  =n/*0 


t 
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x  =  u/V,  ,  v  =v/V,  and  y  v 

*ov/ 

The  velocity  components  may  be  computed  from  the  stream  function  using  Eqs.  (2) 
and  (3)  as. 


u(x,y)  =  |"J?  (x,y) 
v(x,y)  =  fjjr  (*,y) 


(12) 

(13) 


The  system  of  equations  is  closed  with 


dVdx2 


V2* 


3/2  “12^  O 


[u2(x,q(x))  +  v2(x,ri(x))] 


(14) 


1  V>o 

Notice  that  the  only  parameter  is  the  Weber  number,  We  =  g  P/  g  ,  which 

appears  in  the  boundary  condition,  Eq.  (14), 

In  rder  to  solve  the  problem  analytically  we  would  have  to  get  a  solution  of 
the  Laplace  equation,  Eq.  (10)  that  satisfies  the  boundary  conditions,  eq.  (11), 
where  the  free  surface  depression,  T](x),  comes  from  Eq.  (14).  Notice  that  the 
solution  of  Eq.  (14)  is  linked  to  the  solution  of  \|f  through  u  and  v,  thus  the  problem 
is  highly  nonlinear. 

Fortunately  we  use  an  asymptotic  expansion  technique.  First,  we  may 
expand  all  the  dependent  variables  in  terms  of  the  We  number  obtaining: 


n 


q(x)  =  I  We*  ^(x)  +  CKWen+1) 
i=0 


(15a) 
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n 


y(x,y)  *  £  We*  y^x.y)  +  CKWen+1) 
i=0 

n 

(15b) 

u(x,y)  =  £  We*  Uj(x,y)  +  0(Wen+1) 
i=0 

n 

(15c) 

v(x,y)  =  £  Weivi(x)  +  (XWen+1) 
i=0 

Next  we  substitute  Eq.  (15b)  into  Eq.  (10)  to  obtain, 

(15d) 

V2y  =  vVQ  +  WeV2y1  +  We2V2y2  + ...  =  0 

(16) 

Equation  (16)  holds  for  all  We,  and  thus  all  must  independently  satisfy  Laplace 

equations: 

o 

it 

< 

(17a) 

o 

II 

> 

(17b) 

o 

II 

> 

(17c) 

Using  the  same  reasoning,  it  is  easy  to  show  that  all  the  homogeneous  boundary 
conditions,  Eqs.  (11a)  -  (11c)  lead  to: 


(x=0,y)  *  0 

iSO 

(18) 

yi(x,y-*«»)  =  0 

iSO 

(19) 

yj(x->oo,y)  *  0 

iSO 

(20) 

Let  us  next  expand  the  boundary  condition  in  Eq.  (lid): 


(21) 


*¥0 

dx 


d¥i  o  ^¥2 

(x,y=0)  +  We  -g^r  (x,y=0)  +  We2  -g*  (x,y=0)  + 


-1 


Obviously  the  zeroth  order  term  must  be  equal  to  -1  and  all  other  terms  must  be 


zero.  Thus, 


a^(*oM»  = 


-1 


(22) 


and 


-^i(x,y=0)  =  0  i£l  (23) 

The  boundary  condition  in  Eq.  (lie),  and  the  interfacial  jump  contidion,  Eq.  (14), 
require  special  attention.  Using  Eq.  (15b),  Eq.  (lie)  can  be  rewritten  as: 

0  *  y(x,Ti  (x))  =  y0  (x,t|(x))  +  We  y1(x,Ti(x))  +  WeV2(x,ri(x))  +  ...  (24) 


Performing  an  expansion  in  terms  of  We  in  the  neighborhood  of  q  =  0,  we  obtain: 

ety  1 

0  =  vo  (x,0)  +  -=£  (x,0)tj(x)  +  2~^  (x,0)n2(x)  +  ... 


I  av,  1 dV  2 

+  Wei Vi  (x,0)  +  -  jy  -  (x,0)  tj(x)  +  2 (x,0)n  (x)  +  ... 


4*  ... 


(25) 


Substituting  Eqs.  (15a)  and  (12)  into  Eq.(25),  and  rearranging  (ie,  collecting  terms 
We1),  we  obtain: 


0  *  v0  GM»  +  [uo(x,0)  +  ¥1  (*>°)]  We 


+  We^u^  + 1  +  u2ti1  +  u^g  +  ^ri^  + 1 rii  +  V3<x' 

Thus,  the  boundary  conditions  for  x  >  1  are 


if  1  ^1  2  13 

|Vl2  +  2  ‘^■^1  +  “2^1  +  Vb  +^-^2  +  6^1 


,0)j+ We4 ...( 


Vo(x,0)  =  0 


V^x.0)  = 


f  1  2' 

y2(x,0)  =  -  ju^  +  u0il2  +  2Hf^l 


f  1  2  ^0  1  ®  uo  3 

y3(x,o>  =  -  (u0ti3  +  Uln2  +  +  2-^-ni + -dfW  +  6~djr^i 


(27a) 

(27b) 

(27c) 

(27d) 


Notice  that  has  homogeneous  boundary  conditions  for  i  ^  1  except  for  x  ^  1 

where  its  value  is  given  by  Eqs.  (27).  This  set  of  equations  is  the  result  of  the 
asymptotic  expansion  of  Eq.  (lie).  In  order  to  get  a  closed  set  of  equations,  let  us 
focus  our  attention  on  Eq.  (14),  the  interfacial  jump  condition.  We  first  expand  the 
left  hand  side  of  Eq.  (14)  and  then  the  right  hand  side. 

The  left  hand  side  of  Eq.  (14)  can  be  rewritten  as: 

j2  f  n 


~2  X  T|j(x)  We*  +  0(Wen+1) 


d\ 

dn2  2  3 

tir  .  _ *  itr.A  .  Ttr.O 

3  “^1 

^11 

f  d  n3 1 

-  2  +  j  2  "e  +  »  2  ”e  +vve 

dx  dx  dx 

'2d,2 

[dxj 

+  d*sJ 

d  u1dn1dri2  3  d^T] 

l2 

dx2  dx  dx  ’ 2  dx1 

2 

dx  I 

+  0(We  ) 
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Similarly,  the  right  hand  side  of  Eq.  (14)  can  be  rewritten  as, 
We[u2(x,Tj(x))  +  v2(x,tj(x))]  = 

f  du 

We '  uo(x,0)  +  -jg*  (x,0)  tj(x)  + 1  (x,0)  q2(x)  +  . . . 


du.  i  dV  „ 

+We(u1(x,0)  +-^-(x,0)  q(x)  +  |-^<x,0)n2(x)  +  ...)J 


jvo(x,0)+ 3^r(x,i 


0)  T|(x)  + 1^|  (x,0)  ti2(x)  +  ... 


(  ^  ^ 

+  We|vj(x,0)  +  (x,0)n(x)  + 1  (x,0)  q2(x)  +  . . . 


-2 


A* 


(29) 


Using  Eq.  (15a)  we  obtain: 

We[u2(x,q(x))  +  v2(x,t](x))] 

ff  du 

1  {[u.(l<-0>  +  (%  +  Weri1  +  We\  + ...)  + ... 

\i2 


=  We 


+  We 


Ui(x,0)  +  (x,0)  (q0  +  Well  x  +  We2ri2  +  ...) 


dv 


vo(x,0)  +  0^  (x,0)  (ti0(x)  +  Wer^  +  We2Ti2  +  ...)  +. 


+  We 


dv 


vl<x.°)  +  (x»0)  (q0  +  Wet!  x  +  We2^  +  ...)  +  .. 


\ 

2 

+  ... 

;  -1 

(29b) 


Finally,  we  obtain 


iiiBPBSPPIwp 


U 


(u2  +  V2)  We  =  (u20  +v*)  We  +  We2^^ *li  +  ^  +2v0  [d^  vi +  vi  ] 


+  We' 


3 


P*0  f  0  k 

[a^i+vij  +2yo[if^i+df 


l32Vo  2 

9y  T12+2^2  ^l  +  v2 


l]+[^vl+ul] 


i  a2u0 

+  2^0 [if^l + -ft  *2  +  +  "2JJ 


dui  au0 


LK  1 

3vi  dv0  1  d2vQ  2  2 

2['5T’1i  +  vi 

ay  ^l+^y  ^2  +  2^2  ^1  +  V 

+  2vo[  dy  T12+2^y2  ^1+  0y  T11+  ay  T13+  ^2  ^l^  +  e  dy  ~7U  +  V3j 


avj  1a2v1  2  av2  dv0  a2v0 


ia3vo  3 


+  2u„ 


+  2 


dUi 


l^i  2  3*2  ^  Uo  l^O 


dy  ^2  +  2^2  T1l+  dy  ’U  +  dy  ^3  +  "^2^2  +  6  ^3  ^l  +  ^J 


dui  au0  l.  a2uQ  2 

■5r1'i+ui 

dy  ’11  +  3y  ^2  +  2  0y2  Tll  +  ^j 

+  ... 


(30) 


Comparing  the  expression  for  the  left  hand  side  of  Eq.  (28)  with  the  right  hand 
side  of  Eq.  (30)  we  obtain: 


A 

dx2 


=  0 


d2^l  ,  2  2, 

^2-  =  -(U0  +  v0) 


d\ 

dx2 


'du  > 

(3v0  Y| 

[2“. 

■^T^i  +  Ui 

V  > 

+  2v0 

[V^JJ 

(31a) 


(31b) 


(31c) 


The  boundary  conditions  for  these  ordinary  differential  equations  are: 
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T|j  (x-**>)  =  0  (32a) 

dtij 

(x  >o°)  =  0  (32b) 

Equations  (27)  and  (31)  are  the  expanded  forms  of  Eqs.  (lie)  and  (14),  respectively. 
Equations  (17),  (18),  (19),  (20),  (22)  and  (23)  complete  the  set  of  equations  to  be 
solved. 

ORDER  ZERO  SOLUTION 

In  the  previous  section  we  established  a  set  of  infinite  coupled  equations 
for  q  j(x).  In  this  section  we  will  analytically  solve  the  order  zero  equation  (i.e., 

the  functions  qo,  \yo,  uQ,  vQ).  The  equation  for  qo  is  Eq.  (31a), 


with  the  boundary  conditions 


qo(x-*»)  =  0 

(33a) 

d\ 

*f<~->  =  0 

(33b) 

The  solution  is 

o 

ii 

/ — . 

X 

N - ' 

o 

fT 

(34) 

The  differential  equation  for  \|/0  is  Eq.  (17a), 

vV0  =  o 


and  the  corresponding  boundary  conditions,  Eqs.  (18),  (19),  (20),  (22)  and  (23)  are: 
Y0(x-*»,y)  =  0  (35a) 


'.^,«»«-*i:^J!v-  ';  rJ-W  -4 * 9|! 
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Y0(x,y-*~)  =  0 

(35b) 

(x=0,y)  =  0 

(35c) 

d\|/ 

aT^y=<»=-i 

0  <x<  1 

(35d) 

d\lfn 

-§j-(xor=0)»0 

x  >  1 

(35e) 

It  is  convenient  to  use  complex  variables  to  solve  the  two-dimensional  problem. 
Let  O  be  a  nondimensional  analytic  function  given  by, 

O(z)  =  <Kx,y)  +  i  \j»(x,y)  (36) 

and 


z  =  x  +  iy 


and  <J>,  V  are  the  real  and  imaginary  parts  of  O(z)  respectively.  The  function  0(z) 
is  the  complex  potential,  0  is  the  potential  function  and  \ f  is  the  steam  function. 
The  velocity  components  of  u,  v  are  given  by 


u“dx~dz 


u 

dy  dx 


(37a) 

(37b) 


and. 


dO 

dz  SSU'1V 


(38) 


We  recall  the  boundary  conditions  (35d&e), 


u 


0  <x<  1 
x>  1 


A  function  of  x  and  y  that  verifies  these  boundary  conditions  is: 

This  expression  is  the  imaginary  part  of  the  following  analytic  function 
“  [log  (z-1)  -  log  (z+1)] 


(39) 


(40) 


(41) 


We  may  use  Eq.  (38)  to  calculate  the  velocity  components 

u0  =  Re(0)  (log  [y2  +  (x-1)2]  -  log  [y2  +  (x+1)2])  (42) 

v0  =  -Im(<S)  =  i(atg^-atg^j  (43) 

Table-I  shows  the  velocity  of  the  zeroth  order  solution  as  well  as  the  derivatives 
needed  for  solving  Eqs.  (27)  and  (31).  We  now  have  the  solution  of  the  problem  to 
order  zero.  In  order  to  solve  the  problem  to  order-n,  we  must: 

(1)  Compute  u0  (x,y=0)  and  vo(x,y=0)  using  Eqs.  (42)  and  (43). 

d^i 

(2)  Solve  jj-  as  a  function  of  uo(x,0),  vo(x,0)  (Eq.  (36))  with  the  the  proper 
initial  conditions  (Eq.  (32)). 

(3)  Compute  v1(x,0)  from  Eq.  (27b)  for  x  >  1  using  tj1(x),  uo(x,0),  vo(x,0). 

(4)  Solve  V2^  j  »  0  using  y1(x,0). 

<fyi 

(5)  Compute  ux(x,0)  =  (x,0),  v^x.O),  =  -  (x,0) . 
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(6) 


Evaluate 


dx2 


using  Eq.  (31c). 


(7)  Compute  V2(x»°)  fr°m  Eq.  (27c)  for  x  >  1,  using  t\1,  q2,  u0*  v^,  ux,  vv 

(8)  Solve  V2y2  =  0  using  y2(x,0). 


RESULTS 

Figure-3  shows  the  analytical  results  for  the  solution  of  zeroth  order 
transverse  velocity,  u*  computed  at  y  =  0. 

We  have  shown  that  T]0(x)  is  zero.  The  results  for  ,n1(x)  are  shown  in 

Figure-4.  This  is  the  first  correction  term  in  the  shape  of  the  interface.  Using 
q1(x)  and  uo(x,0)  we  can  compute  v1(x,0)  from  Eq.  (27b),  the  non-homogeneous 
boundary  condition  for  \|fj(x,y).  Figure-5  shows  the  curves  of  constant  obtained 

with  PHOENICS  (Spalding,  1991].  Now,  we  have  the  problem  solved  up  to  first 

order.  We  could  compute  the  position  of  the  surface  in  two  possible  different  ways. 
First,  from  q^x),  and  secondly  from  Eq.  (15b) 

V0(x,y)  +  We  V^x.y)  =  0  (44) 


It  is  interesting  to  note  that  these  do  not  produce  the  same  result.  The  reason 
being  that  we  did  not  solve  the  complete  system  of  equations.  In  other  words,  if  we 
compute  the  shape  of  the  interface  as, 

oo 

q(x)  =  X  WeHi^x)  (45) 

i=0 


where  q(x)  will  be  equal  to  the  y(x)  obtained  from  the  equation 
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£  WeVj(x,y)  =  0  (46) 

i=0 

However,  if  we  truncate  the  series  and  we  compute  only  a  finite  number  of  terms, 

we  will  not  get  the  same  result  for  the  surface  from  the  two  methods.  From  the 
result  of  xj/jte.y)  we  can  compute  the  velocity  components  u1(x,0),  v1(x,0).  Figure-6 

shows  u^x.0)  and  v1(x,0). 

Figures  7-8  show  the  contour  maps  for  y2  and  \y3,  respectively.  Figure-9 
shows  the  velocity  components  corresponding  to  the  stream  function  y2.  Figures 
10  and  11  show  the  second  order  eigenfunctions  in  T|,  i.e.,  t)2  and  q3.  Figure-12 

gives  q  x),  the  shape  of  the  interface,  for  five  different  Weber  numbers  (ie,  We  = 
0.1,  0.2,  O.3.,  0.4  and  0.5).  We  arbitarily  define  the  width  of  the  gap  between  the 
interaces  (8)  as  the  width  corresponding  to  half  the  depth  of  the  depression. 
Figure-13  defines  the  gas  gap  width  (8).  Figure  14  depicts  the  width  of  the  gas  gap 
as  a  function  of  the  Weber  number  obtained  from  Figure-12  using  the  definition 
from  Figure-13. 

SUMMARY  AND  CONCLUSIONS 

This  paper  presents  the  results  of  an  analysis  based  on  asymptotic 
expansions  for  the  surface  depression  produced  by  a  plunging  liquid  jet.  The  first 
three  terms  of  a  Taylor  series  expansion  in  the  Weber  number  have  been  obtained 
for  the  surface  depression.  This  approximation  of  the  surface  depression  gives 
correct  values  for  small  and  moderate  Weber  numbers  (ie,  We  <  1).  However,  for 
We  greater  than  unity,  higher  order  terms  become  important  and  the  analysis 
presented  here  is  no  longer  valid. 

The  results  described  in  this  paper  show  that  as  the  Weber  number  is 
increased,  the  terms  that  gain  importance  (i.e.,  the  terms  of  higher  order) 
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correspond  to  a  surface  depression  that  is  increasingly  narrower  in  the  horizontal 
direction.  For  a  Weber  number  of  the  order  of  unity  the  surface  tension  is  no 
longer  strong  enough  to  keep  the  system  stable  and  an  instability  leads  to  air 
entrainment.  That  is,  for  values  of  the  surface  tension  going  to  infinity  (i.e., 
Weber  number  going  to  zero),  the  slope  of  the  surface  is  very  small,  however,  as 
the  Weber  number  is  increased  the  slope  also  increases.  For  a  critical  value  of  the 
Weber  number,  the  slope  of  the  surface  is  such  that  the  surface  tension  is  not 
large  enough  to  keep  the  pool  surface  from  touching  the  plunging  liquid  jet  and 
thus  air  entrainment  is  produced. 

The  results  presented  in  this  paper  cannot  predict  the  route  to  instability 
that  produces  air  entrainment  for  two  reasons.  First,  it  is  not  clear  that  the 
approximation  of  the  surface  depression  is  valid  for  Weber  numbers  that  are  high 
enough  to  produce  air  entrainment.  Second,  the  shape  of  the  surface,  T|(x),  has  to 
be  a  monovalued  function  of  x  and  it  appears  that  at  the  point  where  air 
entrainment  is  produced,  for  every  x  (horizontal  position)  there  is  more  than  one 
t\  (vertical  position  of  the  surface  depression). 

It  appears  that  it  would  be  useful  to  compute  the  surface  position  using  an 
appropriate  multidimensional  Computational  Fluid  Dynamics  (CFD)  tool  having 
surface  tracking  capability.  This  will  allow  the  relaxation  of  the  small  and 
moderate  Weber  numbers  assumption. 
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Figure  1 


Deformation  of  the  pool  surface  produced  by  a  plunging  liquid  jet 
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Liquid  Jet 


Figure  S  Contour  plot  of  the  stream  function,  Vj(x,y) 


Liquid  Jet 


Figure  7  Contour  plot  of  the  stream  function,  \y2(x»y) 


Liquid  Jet 


Figure  8  Contour  plot  of  the  stream  function,  y3(x,y) 
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Shape  of  the  interface  for  different  Weber  numbers 


Figure  13  Definition  of  the  gas  gap  width,  5 


We 
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width  5  as  a  function  of  the  We  number 


